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ABSTRACT 

I present results of simulations of the q=10 and q=20 2-d Potts models in the transi- 
tion region. The asymptotic finite size behavior sets in only for extremely large lattices. 
We learn from this simulation that finite size scaling cannot be used to decide that a 
transition is first order. 
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1. Introduction 

This talk addresses the question of deciding whether a transition is first order 
or not, using Monte Carlo simulations. 

• Very strong transitions are easy to detect. The system behavior is very close 
to the infinite volume behavior. Ergodicity is broken. Thcrmodynamical quantities 
are discontinous at the transition point, with metastable branches. A starting 
configuration half ordered, half disordered will relax to very different states on both 
sides of the transition. 

With extreme statistics, one is able to sample the coexisting phases (although 
this may be forbiddingly costly). The time evolution of any thermodynamical quan- 
tity then shows flip-flops between the phases, and the corresponding probability 
distribution is made of well separated peaks. 

• In less clear cases, one must simulate systems of increasing volumes L d and 
try to convince oneself that the above described very large volume behavior is 
approached. Let me use the language of energy driven transitions in what follows, 
and introduce the energy probability distribution Pl(E). In the transition region, 
it has two peaks of heights P£ and Pf, separated by a minimum of height p™ m . 
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For a first order transition, at fixed P£/P^, one has 



^ pmin 

— r In — ^ > L -> oo 

This has been proposed by Lee and Kosterlitz 1 as an indicator of first order phase 
transitions (see also 2 3). 

Another class of indicators are moments of the distribution Pl(E), that goes to 
zero in the large volume limit, for all temperatures, but at a first order transition 
point. An example is the energy fluctuation CV/L d = (3 2 (< E 2 > - < E > 2 ), 
another is Binder's famous cumulent BL = 4(1— < E A > / < E 2 > 2 ). One can 
also plainly look at a plot of Pl {E) and decide "by eyes" whether it approaches two 
delta functions as L grows. 

• The most sophisticated (and trustworthy) method is finite size scaling. One 
insist in seeing the finite size behavior as predicted by the theory in the vicinity of a 
first order phase transition. In that case only can one be pretty sure that the trend 
observed for lattices of increasing sizes does continue up to the thermodynamical 
limit. One insist on seeing 

^ln P , { V ~A/L 
L d P°(L) ' 

CV max /L d + CV^/L d 

where CV max is the maximum of the specific heat, 

BL mm /L d ~ BL« + BL^/L d 

where BL m i n is the minimum of BL 

My interest in first order phase transitions started with the 1988-1989 contro- 
versy over the order of the deconfinement phase transition of pure SU(3) gauge 
theory, and the question whether some modified 3-d Z(3) Potts model has a first 
order transition or not. Our data 4 were very convincing showing that BL min has 
a non zero large volume limit, however we did not observe the predicted l/L 3 finite 
size behavior (our data behave nearly like l/L 2 ). This led us to investigate the 2-d 
q=10 Potts model, as an example of a model with a strong first order transition, 
searching for the predicted finite size behavior. 

2. Exact Results 

Those have been obtained 5 for models that can be represented by a contour 
expansion with small activities, like 6 the q states Potts model for large q. In such 
a case, the partition function for a L d lattice with periodic boundary conditions can 
be written as 



Z(0, L) - e - Ld ^ + qe - Ld Pf°W + 0{ t,-hL )e - &m L* ■ b > (2.1) 



where / (/3) and fd{0) are smooth L independent functions. The free energy is 
f(P) = min{f (/3), fd(/3)}. The phenomenological two gaussian peak model of the 
energy probability distribution Pl(E) introduced by K. Binder and D. Landau 7,8 
follows through inverse Laplace transform. The above exact result fixes the relative 
weights of the two peaks: At the infinite volume limit transition point, [3 = (3 t , the 
ordered and disordered peak weights are exactly in the ratio q to one. 

The two gaussian peak model is not however a good representation of Pl (E) for 
all E's. It fails to describe the region between the two peaks, and does not account 
for the observed L dependence of the position of the two maxima of Pl (E) when 
P = (3f This would require the understanding of the correction term in Eq.( 2.1.) 

To the order in 1/L d we consider, all quantities are expressed in terms of /3 t 
and of the energies and specific heats of the two coexisting phases. The transition 
temperature, E Q , Ed and the difference C — Cd are known exactly for the 2-d Potts 
models 9. It follows fromEq.( 2.1) that the specific heat 

CV = (3 2 L d (< E 2 > - < E > 2 ) 

has a maximum at 

P(CV max ) = (3 t - + % + Od/L^). 

The height of this maximum increases linearly with L d 

CV max = L d ®-(E - E d ) 2 + CV^ + (D(l/L d ). 

whereas for fixed (i ^ (i u CV{j3) goes to a constant, as L goes to infinity. One finds 
that BL reaches a minimum equal to 10,5,11 
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at the point 
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Expressions of the coefficients (^\, BL^ 2 \ fl^y and CV^ as functions of the 
Ei's and d's can be found in 11. Although the use of BL m i n as an indicator of 
the order of phase transitions has been much publicized, CV max /L d is as good an 
indicator indeed. The value of BL depends on the choice made of the arbitrary 
constant one can add to the definition of the energy, this leads to introduce the 
quantity 12 

<(E-<E>) 4 > 



which is independent of such a constant. U4 is strictly larger than one, but at a first 
order transition point, in the infinite volume limit. For large but finite volumes, Z74 
reaches a minimum 



at the point 
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BL min (or CV m ax/L d ) and UA min are dual since BL min (or CV max /L d ) going to 
zero means second (or higher order), whereas UA m i n going to one means first order. 

The above formulae for the extrema CV max , BL min , f74 mi „ and the correspond- 
ing effective /3's have higher power law corrections that may hide the asymptotic 
behavior on lattices that can be simulated. In contrast, the expressions for bulk 
averages evaluated at the (infinite volume limit) transition point (3 = [3 t do not have 
power law corrections, as a consequence of Eq.2.1. The average energy is given by 

and the value of the specific heat is 

CV (P t ) = + j^(Eo E d fft + 0(e-" L ) 

The energy at 3 t does not depend on the lattice size, up to exponentially small 
corrections. This provides an efficient estimator of the transition temperature 6,13 
by the following "two-lattice method" . One simulates lattices of increasing sizes 
Li < L2 < L3 < .., and consider /3 e /j(Lj, Li+i), the solution of the fixed point 
equation 

E Li (0) = E Li+ M 

for i = 1, 2, . . .. The estimate /3 e ff(Li, L i+ i) converges towards (3 t with exponential 
pace. The spectral density (a.k.a Ferrenberg-Swendsen FS, or reweighting) method 
is invaluable for locating extrema and zeros with Monte Carlo data as input. It 
allows to reconstruct the value of a thermodynamical average for any (3 from one 
run performed at a given (3m.c. in the vicinity of the transition. It is well known that 
this method has problems close to a second order point. With moderate statistics, 
it predicts extraneous extreme for e.g. the specific heat 15. This never occurred to 
us with first order points. 

To summarize, the extrema of CV max /L d , BL m i n , and U4 m i n behave in the 
large volume limit like X 1 + X 2 / L d + 0(1/ L 2d ), CV{f3 t ) behaves like X 1 +X 2 /L d + 



0(e bL ). The four different constants {X 1 } are exactly known. One single unknown 
parameter, e.g. the ordered specific heats C Q , fixes the four X 2, s. 

3. Simulation of the q=10 model 

We have performed a w 800 CRAY X-MP hours simulation 12 of the q = 10 
2-d Potts model, in order to determine how large L has to be in order to see the 
asymptotic regime described in 5. This is a model with a strong, obvious, first order 
transition, with 16 a correlation length ~ 6. The precision of our data is 

better by more than one order of magnitude than in 8. We have compared our 
results for the extrema of CV/L d , BL and U4, and for the value CV((3 t ) with the 
large volume predictions. 

We simulated lattices up to L = 50 (where the autocorrelation times are t$ « 
.9 10 6 and t^s ~ 2.5 10 6 ). For all four quantities, we see deviations from the 
X 1 +X 2 / L d limiting behavior. There is nothing to worry about that, it only means 
that our precision is good. Really disturbing however, is that these corrections do 
not seem to behave simply as function of L, and are definitely not under control. 
The values for the four slopes {X 2 } one would infer from our data give inconsistent 
estimates of C a . Three possible explanations are i) Eq.( 2.1) is only proven in the 
large q limit, it may not hold down to q — 10. ii) Much larger lattices may be needed 
in order to extract the true asymptotic behavior, although in our data Pl(E) has 
a textbook first order shape, iii) A programming error is always possible. 

Before doing the simulation, we hoped that CV(j3 t ) would be asymptotic earlier 
than CV max , since corrections are 0(e~ bL ). The data do not substantiate this hope, 
CV((3t) has larger error bars, but does not seem to reach its asymptotic behavior 
earlier. Note that the estimate of C we get from CV(Pt) is much higher than the 
others. 

4. Simulation of the q=20 model using the Multicanonical Algorithm 

The conventional Metropolis (and Swendsen-Wang 17) algorithm suffers from 
exponential slowing down. This makes simulations on lattices much larger than used 
to day, impossible even with vastly more powerful computers. It has been proposed 
by B. Berg and T. Neuhaus 18 to perform the simulation with an Hamiltonian 
designed in such a way that Pl(E) is very smooth between E Q and Ed, and to 
reweight the events when computing expectation values. The new "multicanonical" 
algorithm has only polynomial slowing down. 

We 19 used this algorithm in order to simulate the q=20 Potts model. It has a 
stronger first order transition and a smaller correlation length than the q=10 model. 
This means that the large volume regime sets in for smaller lattices for q = 20 than 
for q = 10. We ran mainly on IBM RS6000 workstations, and simulated lattices as 
large as 38 2 . In contrast with the q=10 case, the values of the four slopes {X 2 } 
one infer from our data give consistent estimates of C . As an example, Figl gives 
our results for CV max /L d together with the theoretical estimate using the value 
C = 5.2 ± .2, and Fig2 gives our results for CV((3 t ) together with the theoretical 



estimate using the same value for C . Note that CV(/3t), reaches its asymptotic 
behavior much earlier than CV max , as predicted by the theory. In conclusion the 
asymptotic behavior predicted by K. Binder, and later proven by C. Borgs and 
R. Kotecky only sets in for very large lattices. The lattice size must fulfil the 
conditions L >> £, L d ~ x >> l/A od where A od is the order-disorder surface tension 
(If Widom's relation 20 holds this condition is equivalent to the first one), and 
L d » C /(E - Ed), L d >> Cd/(E a - Ed), where >> means five to ten times 
larger. It is unfortunate that for such large systems, the transition is blatantly first 
order. 
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Figure 1: CV max /L 2 as a function of 1/L 2 for the 2-d q=20 Potts model. 



Figure 2: BL min /L 2 as a function of 1/L 2 for the 2-d q=20 Potts model. 
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